NONFACTORIZABLE SOFT CONTRIBUTIONS IN THE B > r]cK, XcqK DECAYS WITH THE 

LIGHT-CONE SUM RULES APPROACH 

Zhi-Gang Wangi'2,3 *^ ^nd Tao Huangi-^ 

1 CCAST (World Laboratory), P.O.Box 8730, Beijing 100080, P. R. China 
2 Institute of High Energy Physics, P.O.Box 918 ,Beijing 100039,P. R. China 
^ Department of Physics, North China Electric Power University, Baoding 071003, P. R. China 

In this article, we calculate the nonfactorizable soft contributions in the B VcK, XcoK 
' decays with the light-cone QCD sum rules approach. Our results show that the nonfactorizable 

I corrections from the soft gluon exchanges in the decay B rjcK are of (15 — 30)% and should be 

^vQ ■ taken into account. As for the decay B — » x^oK, the factorizable contributions are zero and the 

nonfactorizable contributions from the soft hadronic matrix elements are too small to accommodate 
the experimental data. 
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> 

Q> . I. INTRODUCTION 

(N 

, The nonleptonic decays of the B meson have attracted much attention in studying the nonperturbative dynamics 

o- 

, of QCD, final state interactions and CP violation. The exclusive B to charmonia decays are important since those 
Oh! 

' , decays B J/ipK^ricK^XcjK are regarded as the golden channels for the study of CP violation. The quantitatively 
M-i. 
D . 

(— I , understanding of those decays depends on our knowledge about the nonperturbative hadronic matrix elements of the 



X 



operators entering the effective weak Hamiltonian. In Ref. [1], the authors propose an original approach called QCD- 
improved factorization to deal with the two-body nonleptonic decays of the B meson. In this approach, the decay 
amplitudes are expressed in terms of the semileptonic form factors, hadronic light-cone distribution amplitudes and 
hard-scattering amplitudes. The semileptonic form factors, the light-cone distribution amplitudes are taken as input 
parameters and the hard-scattering amplitudes including nonfactorizable corrections due to the exchanges of hard 
gluons are calculated by perturbative QCD. For the exclusive B to charmonia decays , the QCD-improved factorization 
approach is broken down due to the divergence arising from the soft-gluon exchanges, moreover, the theoretical 
branching fractions are too small to accommodate the experimental data [1~5]. In Refs. [2,3], the authors observe 
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that for the exclusive B rjcK decay, the nonfactorizable corrections to the naive factorization are infrared safe 
at leading-twist order, the spectator interactions arising from the kaon twist-3 effects are formally power-suppressed 
but chirally and logarithmically enhanced; for the B XcoK decay, there are infrared divergences arising from 
the nonfactorizable vertex corrections as well as logarithmic divergences due to the spectator interactions even at 
leading- twist order. 

The effects of soft gluons which break down factorization are supposed of order 0{AQCD/'mb) and neglected in the 
QCD-improved factorization studies, however, no theoretical work has ever proved that they are small quantities. For 
the color-suppressed B to charmonia decays, there may be significant impacts of the nonfactorizable soft contributions. 

On the other hand, the QCD light-cone sum rules (LCSR) approach provides a powerful tool for calculating the 
exclusive soft hadronic amplitudes [6-8]. The LCSR approach carries out the operator product expansion near the 
light-cone a;^ ~ instead of the short distance a; « while the nonperturbative matrix elements are parameterized by 
the light-cone distribution amplitudes which classified according to their twists instead of the vacuum condensates. 
For detailed discussions of this approach, one can consult Ref. [9] . The LCSR approach has been applied to study the 
nonfactorizable hadronic matrix elements due to the soft gluons exchanges and gave satisfactory results [10-12]. 

Experimentally, the B rjcK decay was observed by CLEO, BaBar, Belle Collaborations and the B XcoK 
decay by Belle Collaboration with relatively large branching fractions [13,14]. The large discrepancies between the 
theoretical and experimental values for those decays call for considerations of new ingredients and mechanisms. 
It is interesting to study the nonfactorizable soft contributions in the pseudoscalar and scalar charmonia decays 
B — > rjcK, XcoK with the LCSR approach. 

The article is organized as follows: the factorizable contributions from the effective weak Hamiltonian are derived 
in Sec. II; the soft hadronic matrix elements {'qcK\0\B) and {xc'aK\0\B) are calculated with the light-cone sum rules 
approach in Sec. Ill; numerical results are presented in Sec.IV; the section V is reserved for conclusion. 



The effective weak Hamiltonian for the h — > sec decay modes can be written as (for detailed discussion of the 
effective weak Hamiltonian, one can consult Ref. [15]) 



II. EFFECTIVE WEAK HAMILTONIAN AND FACTORIZABLE CONTRIBUTIONS 




(1) 
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where V^j's are the CKM matrix elements, Cj's are the Wilson coefficients calculated at the renormalization scale 



^ 0{mh) and the relevant operators Oi are given by 



Oi = {Sabi3)v-A{Cl3Ca)v-A, C>2 = {Saba)v-A{Cl3C0)v-A, 

C'3(5) = {Saba)v-A'^{q(3q(})v-A{V+A), ^'4(6) = {Sabi3)v-A^^{q0a)v-A{V+A), 



o. 



7(9) 



-{saba)v-A^^eq{qf^qi3)v+A{v-A), 08(10) = :^(sa^/3)y-A^eg(g^ga) 



y+A(y-A)- 



(2) 



We can reorganize the color-mismatched quark fields into color singlet states by Fierz transformation, (for example, 
Oi = 7^C2 + 202, Nc is the color number and taken as 3.) and express the effective weak Hamiltonian Hyj in the 
following form. 



-VtbVt*s 



C3(/X) + O3 + 2C4(m)03 



C9(m) + ^^)O9 + 2Ci0(m)O9 



(3) 



where 



^2(3,9) = (C7/.(1 - 75)c)(s7^(l - 75)&) , ^2(3,9) = (c7^(l - 75 ) y c) (57" ( 1 - 75)^6) , 
©5(7) = (C7/.(1 + 75)c)(s7'*(l - 75)&) , ^5(7) = (57^1 + 75) yc)(s7''(1 - 75) Y&) , 



(4) 



here A"'s are SU{3) Gell-Mann matrices. 

The factorizable matrix elements of the operator Oi for the decay B — > ricK can be parameterized as 



{Vc{p)K{q)\HMP + 'l)) = ^ 



3 J "V 3 

(^cb)|c7^(l - 75)c|0)(if («)|S7'^(1 - 75)6|B(p + q)) - Vt^VZ (C^ii^) + Cr{fx) 



+ 



{Vcip)\cjf.{l + ^5)c\0){Kiq)\sj^{l - j5)b\Bip + q)) 



+ 



C4(M) - - g8(M) + Cio(M) -| I (,,(p)|e7,75c|0)(i^(.)|^7^6|E(p + ,)) , (5) 



3 



where the meson momenta are explicitly specified and chosen as = m^^ . The r^c meson decay constant is defined 
by the relation, 

(r?c(p)|c(0)7^75c(0)|0) = -i/^^p^. (6) 

The decay constant /^^ can be estimated from the QCD sum rules approach with the current = c'yalbC [16] 
or phenomenological potential models, in fact, the values obtained from those approaches do not differ from each 
other much; new estimation based on the nonperturbative approach of coupled Schwingcr- Dyson equation and Bethe- 
Salpeter equation is in preparation. We use the value obtained from the potential model in this article [17]. 
The B — K form factor can be parameterized as 

{K{q)\-s^,h\B{p + q)) = {2q+p),F+^{p^) - !!^^^<p^(F+^ip^) - F°^{p'')), (7) 

the above form factors Fgj^{p'^), F'^^[p^) can be estimated from the light-cone sum rules approach [18-20], here we 
take the value 

^^bkKc)= 0.42 ±0.06. (8) 
The concise expression for the factorizable matrix elements in the decay B ijcK can be written as 
{V,{p)K{q)\H^\B{p + q)) = ^ I^V^.V;, (^C2(m) + ^) 

There are no factorizable contributions in the decay B XcoK as the vector current C7^c is conserved and has a 
vanishing matrix element with the XcO meson, 

(Xco(p)|cc|0) = /x,o"^XcO > 

{Xco{p)K{q)\H^\B{p + q)) = 0, (10) 
here /p^^g and m^^^ are the decay constant and mass of the XcO meson respectively. 

III. LIGHT-CONE SUM RULES FOR {ricK\d\B) AND {xcoK\d\B) 

In the following, we apply the approach developed in Ref. [10] for the B ^ nn channel to estimate the con- 
tributions from the soft-gluon exchanges in the B r]cK,XcoK decays. Firstly, let us write down the correlation 
functions, 
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q, k) = e j d^a^e-^ff-")- j d^e'^^-^'^y {<:)\T{f^^{y)0{Q)ji {x)}\K{q)) , (11) 
F^^,{p,q,k) = *2 I d^xe-'^^-i^^ j d^ye^^P-^'^y{Q\T{j^'\y)d{Q)ji{x)}\K{q)) , (12) 

where j^" = 07^750, j^"" = cc and = mbbi^^u are currents interpolating the rjc ,XcO and B meson fields, respectively. 
We take the same notation O for the operators Oj as our final result indicates that they have the same analytical 
expressions due to their special Dirac structures for the decays B — > VcK, XcoK . 

The correlation functions can be calculated by the operator product expansion method near the light-cone x'^ ~ 
j/^ ~ {x — y)"^ ~ in QCD. They are functions of three independent momenta chosen to be g, p — fc and k by 
convenience. Here we introduce the unphysical momentum k in order to avoid that the B meson has the same four- 
momentum before {p — q) and after the decay [P) . In such a way, we can avoid a continuum of light contributions in 
the dispersion relation in the U-channel. The independent kinematical invariants can be taken as {p — q)^, {p — k)^, 
q^ , fc^, = {p — k — qY and p^. We set fc^ = and take = m% = 0, neglecting the small corrections of the 
order 0{m\ /rn?g). The momentum p^ is kept undefined for the moment in order to make the derivation of the sum 
rules without restriction. Its value is going to be set later in this section, and chosen p^ = ni^^ , rn^^^ . The values of 
(p — fc)^, [p — q)^ and P^ should be spacelike and large in order to stay far away from the hadronic thresholds in the 
B, ric and XcO channels. All together, we have 

q^=k^ = 0, p'^ (undefined), \{p - k)\'^ > Aqcd,\{p - q)f » ^qcd, \Pf » ^qcd ■ 

The decomposition of the correlation function in Eq.(ll) with the independent momenta are straightforward and it 
can be divided into the following Lorentz invariant amplitudes*, 

F^^ {p, q, k) = {p- k)pFr,^ {p, q, k) + q^F^^ + k,F^^ + Cp^x^qf^p^k^F^' . (13) 

According to the basic assumption of current- hadron duality in the QCD sum rules approach [16], we insert a complete 
series of intermediate states satisfying the unitarity principle with the same quantum numbers as the current operators 
ip°) j'^"" and into the correlation functions in Eq.(ll) and Eq.(12) to obtain the hadronic representation. After 
isolating the pole terms of the lowest pseudoscalar rjc and XcO mesons in the charmonium channels, we get the following 
result. 



Fr,^ and F^^^ are functions of q, p — k and k, their momenta dependence will be written explicitly when necessary, for example, 
FriciP ~ ^) FxcoiP ~ k), with emphasis on the dependence on the momentum p — k. 
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where 



n,.(xco)((f^-a)',i^',p') =i I d4xe-^(^'-'^)-(77e(Xco)(p- A:)if(-g)|r{O(0)jf (:r)}|0). (16) 

In above equations, p^p" and Sq'' are the spectral density and threshold parameter of the lowest excited resonances and 
continuum states in the r/c channel, respectively; while p^''° and Sg"" are the corresponding ones in the XcO channel. 
In the limit of large spacelike momentum {p — fc)^ <S; m^^ , rri^^^ , the correlation functions in Eq.(ll) and Eq.(12) can 
be calculated in QCD at the level of quark-gluon degrees of freedom and rewritten in the following forms by applying 
dispersion relation. 



1 r ^ lvcisF^p^{F^J{s,{p-q)\P\p^) 



We can approximate the hadronic spectral densities pj^= and p^''^ above the thresholds of the lowest excited resonances 
and continuum states by the corresponding ones from QCD calculations with the assumption of quark-hadron duality, 

p^^(p>^-)(s,(p-g)^p^p2)0(,_,g=(,xco)) = -Im,F;=(^^^„J(s,(p-9)^p^p^)e(s-s3^(sr))• (is) 

Now we explore the analytical properties of the amplitudes n^^((p— g)^, P^,p^), n^^o((p— g)^, P^,p^) in the P-channel 
and insert a complete series of hadronic states with the same quantum numbers as the B meson into the correlation 
functions in Eq.(16). After isolating the lowest pole terms of the B meson contributions, we obtain the following 
results, 

97)2 f poo pB tgl p2 2\ 

n..(xco)((p - g)^ p\v') = _ (^^/^^2 _ ^, (^c(xco)(p - k)K{-q)\omB{p + g)) + j^^ ds' ■ (19) 

The hadronic spectral densities p^^ and p^^^ above the threshold of the continuum Sq can be approximated by the 
corresponding ones at the level of quark-gluon degrees of freedom. 

Finally, we obtain the correlation functions in hadronic representation. 



ifvo{p-k)p fBTTl- 



fn^^ — [p— ky — ie m\ — {p — q)^ — 



-{ricip - k)K{-q)\OmB{p - q)) + ■ ■ ■ , 



IS 



= {p-k)pFr,^{p,q,k) + ---; (20) 



(0|j^ -°(0)IXco(p - fc)) Mi^^/HM - 'Zlil'WiO) 

— (p — — i£ 



- ^.^^ _ _ _ ^1 _ (, _ - fc)i^(-.)P(0)|B(p -,)) + .... (21) 

Here we do not show the contributions from the higher resonances and continuum states above the thresholds explicitly 
, they can be written in terms of dispersion integrals and the spectral densities can be approximated by the quark- 
hadron duahty ansatz. In Eq.(20), we select the relevant terms with tensor structure {p — k)p, which corresponding 
to the contributions from the pseudoscalar mesons, for example, the r]c and B mesons. 

In order to suppress the contributions from the excited and continuum states in the charmonium channels, we 
can perform n-th derivative with respect to the momentum (p — fc)^ in Eqs. (14-15) to obtain n-th moment sum rules 
for the correlation functions Ylr]c{{P ~ '1)^jP^jP^) and II-)(^^q{{p — q)"^ , P"^ , p"^) in hadronic representation, 

i{p - fc),n,J(p - qf, P',p') = —- / ds'-jf—^^lm,F^^ {s, {p - qf, P^p") , 

1 r^l^ , {ml +Qlr+^ r ds 



f„, J Ami (S + Q 

(m2 +Q2)«+i 



(22) 



n^Aip-<i)\p\p'') = —-^ ° ds rin^^n+i Im.f,.„(,s,(p-g)^p^p^), 



XcO ■' XcO 

1 (m2 +Qlr+^ ds' o o 

/ . f „,, ImJm..F,,„(5,a^P^p^), (23) 



"^xco /xeo Arr.^ + - (p - q) 



where the imaginary parts with respect to s and s' are given by 

Im,Im,,j;=(i^^,J(s,s',p2,p2) = i(p- fc),/,^(/^^„m;,,j7r2/Bm|(r?e(Xco)(p- fc)i^(-9)|5(0)|B(p- 5)) 

^(s'-m|)^(5-m2^(^^„)) + --- . (24) 

Here Qq is the parameter for QCD sum rules in the charmonium channels, the spectral densities above the thresholds 
can be approximated by the corresponding ones from QCD calculation and not shown explicitly. 

The Borel transformations with respect to {p—q)^ in the B channels in Eqs. (22-23) are straightforward. Comparing 
with Eqs. (14-16), we can obtain the following results^, 

o — 1 e m'-^ 

{p-k) =ifr,Xfxco^xco)fBmB{nc{Xco){p-k)K{-q)\0{Q)\B{p-q))—-^ TTWY^^ir^ ' 



^Here we prefer the notations f?)c(p — k), F^^„(p — k) to Fr,^{p,q, k),F^^g(p,q, k), with emphasis on the dependence on p — fc. 
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here Btrans denotes both the n-th derivative and Borel transformation, is the Borel parameter in the B channel. 

The contributions from the excited and continuum states are not shown explicitly for simplicity. 

Then we can analytically continue from the space-like region <C to the time-like region > 0, and 

choose P^ = m%. Now we carry out the operator product expansion near the light-cone to obtain the representation at 

the level of quark-gluon degrees of freedom for the amplitudes Fr,^ and F^^^ . Firstly, let us write down the propagator 

of a massive quark in the external gluon field in the Fock-Schwinger gauge [21,22], 

, 1 



{0\T{q,{x,)q,ix2)}\0) = iJ (^4 e"''^'^""^^ | " / dv G^,'^ {vx, + (1 - v)x2) (^^ 

^ 



-2 (fc^-m^)^ ^^- - W^"^""' - ^^^"H I ' 

here Gf^'^ is the gluonic field strength, gs denotes the strong coupling constant. 

Substituting the above b and c quark propagators into the correlation functions in Eqs. (11-12), we can obtain 
the hadronic spectral densities at the level of quark-gluon degrees of freedom. The following three particle kaon 
distribution amplitudes are useful in our calculation, 

- twist-3 distribution amplitude 

{Q\s{0)ai^uj5Ga(s{vy)u{x)\K+ {q)) = ifsK [{qaqij,g0i^ - H3%9av) 

- qpq.Qa^)] j Vai^au fi)e-">^^"^+y^"^^ ; (27) 

- twist-4 distribution amplitudes 

{0\s{0)i-y^G^^{vy)u{x)\K+{q))=q^^^^^^^^fK [ Pai.^|| (a,, /i)e-^«(-«^+''™3) 

qx J 

+ (5^a9/3 - g^fsqa) J Va^Mocr, ^)e-*«(-«^+«-«=^) ; (28) 
(0|s(0)7^75G„^(^;y)n(a;)|Jr+(g)) =g^^^i^^^^/K / Dai<A||(ai,M)e-'«(""^+'"'"^^ 

qx J 

+ {9Um - g^qo) I T^aiMc^i, ^)e-'^i--^+y^-^) , (29) 



where 



1 A^C"^ 



Vai = daida2da3,5{l-a\-a2-ai), g^/j = gap - . (30) 



The twist-3 and twist-4 light-cone distribution amphtudes can be parameterized as 



't>3K{oii, IJ) = 360Qia2a3 yl + a(/i)-(7Q3 - 3) + 6(At)(2 - 4aia2 - 803(1 - as)) 
+ c(/Lt)(3aia2 - 2(^3 + 3a^)) , 



^_L(ai,/Lt) = 305^(/i)(Qi - 0:2)03 
(^ll(oi,/x) = 120(5^(/Lt)e(/i)(oi - 02)010203 
(j)±{ai,ii) = 30(5^(/i)o3(l - 03) 

(^||(Oi,/U) = -1205^(/z)oiO2O3 



-+2e(M)(l-2o3) 



-+2e(/.)(l-2o3) 



+ e(/x)(l-3o3) 



(31) 
(32) 
(33) 
(34) 
(35) 



Those parameters in the hght-cone distribution amplitudes can be estimated with the QCD sum rules approach 
[9,23,24]. In practical manipulation, we can neglect the a(yK), 6(/Lt), c(/x) and t{p) dependence and the asymptotic 
forms will be taken. 

After carrying out the operator product expansion near the light-cone, we obtain the following expressions for 



the F^^ and F^^^ , 



■m-bf: 



3K 



["I 



4>i.K{ai,^i) 



dx- 



2a;2(l - x) 



47r^ Jq J — (p — g(l — oi))2 Jq — {p — k — va3q)^x{l — x) 

q-{p-k) \{2 -v)q-k + 2(1 - v)q ■ (p - k) 

,2 , 



"^^■^^ /"^^^ /"p^ '^||("i'M) + <^±(a»,Ai) 



4^2 



dx- 



Fx.oip-k) = 



™2 f /•! 

_rrhJK_ 
47r2 

rubmJsK r 



dv / Po, 



2x\l - x)q ■ {p - k) 
''ml - {p- q{l - oi))2 Jo rnl - {p-k- va3q)'^x{\ - x) 

^ , 2x^(1 -x){2v-?,)q-{p-k) 
dx- 



<j)±{ai,fi) 



* ^ (P ^ 9(1 ^ Q^i))^ Jq "fTT-c ~ {P~ k ~ va3q)'^x{l — x) ' 



47r2 



(psKiai,!!) 



^ dx ^^^^ ~ ^^^^ - 2;)(1 - 2?;))g ■{p-k)q-p 



K 



^ dv J Va^ ^2_^p_ qf^i _ ^^^-^2 — ^2 _ _ ^ _ ya3q)'^x{l - x) 

^ ' ^-^ (o,,^) -20^(0,, m)) ^ [x"^ - x{\ - x){\ - 2v))q ■ {p - k) 

' dx- 



47r2 







dv / Va 



ml — {p — q{l — oi))2 Jq m^ — {p — k — va3q)'^x{l — x) 



(36) 



(37) 



Here we will take a short digression to discuss the technical details. In calculation, we will encounter x-integral 
of the form 



/ 



d'^x^^f{p,q,k,x) 
q ■ X 



in the coordinates representation, which can be formally written as 



/ 



d'^x-^f{p, q, k, x) = A{p, q, k){p - k)p + B{p, q, k)qp + C{p, q, k)kp + epi3x^q'^p^k^D{p, q, k) . 



(38) 
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Multiplying both sides of above equation by qp and taking the Chiral limit = m|- = 0, the expression of the 
relevant quantity A{p, q, k) can be obtained. 

It is easy to perform the Feynman parameter x integral in Eqs. (36-37), we prefer this form in order to facilitate 
the Borel transformation and n-th derivative. For the case of massless quark loops, one can integrate out the variable 
X directly. 

In the following, we write down the dispersion relations for the correlation functions at the level of quark-gluon 
degrees of freedom, 

,N mb f\ , T""^ , r , x^{l-.r) S-2-P'- 

PrjcKP- ^) = "'^ / / / dx- 

47r^ Jo J J^2 ^4^2 Jo i 



Si - (P - S2-ip- 
|/3if'/'3if (a,, [^^(P^ + "^1 - «1 - + (1 - V){S2 - P') 

-mbfK{4>\\{ai,lJ') + 4>±{ai,fj,)) - mbfK^^{oci,n){2v - 3)| 

5{m\ — m^^ai - si(l - ai))S{ml — x{l - x)vazP^ — x{l — x){l — vaz)s2) H ; (39) 

^ , ,^ TOfaTOc , , r"''" , , x"^ - X(l - X)(l - 2v) S2 - 

Stt^ Jq J J„2 ^4^2 Jo si-{jp- qy S2-{p- ky 

{hK(l)3K{oii, n){m\^^ - Si) - mbfK{(i>\\ («», m) - 2^_L(Q!i, m))} 

(5(m^ — m^^^ai — si(l — ai))5{ml — x{l — x)va3P^ — x{l — x){l — va3)s2) H • (40) 

Again, the higher resonances and continuum states contributions are not shown explicitly for simplicity, as they are 
Borel transformation or n-th derivative suppressed. Here we can introduce some notations to simplify the cumbersome 
expressions in Eq.(39) and Eq.(40) respectively, 



1 / , L 4to^ \ 1 / , L 4m? 



^^=o l-\ 1--^ , Xf = - l + 



21 V si ' 21V si 

x{l-x){s-P^) SB-m^^{mU 
In performing the 6 functions integrals in Eq.(39) and Eq.(40), we note that in the space-like region <^ 0, 
the condition ao > etc can be warranted. Performing Borel transformation in the B channel and n-th derivative in 
the r]c, XcO channels, then matching with Eq.(25), finally we obtain the sum rules for the nonfactorizable soft matrix 
elements, 

(??e(p - k)K{-q)\d{0)\B{p - g)) = , , 2 n ds r xdx ['da P dl3[f3KhK{l - a,a- 

^T^^fBfr,^m% J 4^2 J^. Ja^ Ja^ 

|(s _ p2yi _ ^) _ (1 _ Kc(l - Q') + Q(s - P^-ml) j 
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-mbfK{4>\\ {l-a,a- P,(3) + 4>±{l-a,a- /?, /?)) 

2ac 

-mbfK<p±{l- a,a- f3,/3){^ -3) 
^""^^ ^ I s + Ql ] ^ ' ^^^^ 



{Xco{p-k)K{-q)\0{0)\B{p-q))= "^'"^^ f"" [' da r dpifaK^Kil - a,a - (3, (3) 



2 2 

''''' ' - m6/x(0||(l-a,a-/3,/3)-2<^i(l-Q,a-/3,/3))] 



{.-(l-«)(l-2^)} 



n+1 



In above expressions, is chosen to be large space-like squared momentum (|P^| ^ m^) in order to stay far away 
from the hadronic thresholds in the channels of the B and charmonia currents, the values of are small positive 
quantities but not always small enough to be safely neglected, we can perform the following approximation for the /3 
integral, 

/" d0G{s, x,a,(3) = l [ - / ' j d(}G{s, x, a, /?), (44) 
here G is an abbreviation for the integral functions and can be written as 

G{s,x,a,(3) = ^(s,a;, a, /3)(/>3if (1 — a, a — /?,/?) + B[s,x,a,(3)(j)\\(l — a, a — (3,(3) + C{s,x,a, (3)<j)±{l — a, a — (3, (3) 
or 

G{s, X, a, (3) = D{s, x, a, (3)((i3k{^ — a, a — /3, /3) + E{s, x, a,P)4>^\ {1 — a,a — P, (i) + F{s,x,a, (i)(p±{l — a,a — P, (i), 

A, B, C, D, E, F are formal notations. We can expand the light-cone distribution amplitudes ^3k, ^||, (t>\\ and (p± 
in terms of Taylor series of /? **, for example, 

d 

03if (1 - a,a- P,(3) = 03k(1 -a, a- P,(3)\fs=o + -g^hK{l - a, a- (3, (3)\i3^o(3 

+\^cl)3K{l-a,a-(3,(3)\p=o(3^ + --- , (45) 



**For very small ttc, we can approximate the integral j^'^ d(3G{s,x,a,(3) by acG{s,x,a,(3)\i3^o, then analytically continue 
into the timelike region, P^ = m%. 
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and analytically continue into the timelike region, = m%, then complete the integral J^" df3G{s, x, a, (3). The 
explicit expressions for the physical matrix elements {r]c{p)K {—q)\O{0)\B{p — q)) and {Xco{p)K{—q)\0{0)\B{p — q)) 
are lengthy due to the re-summation of all the Taylor series of /3, here we show only the leading terms explicitly, 



{il,{p)K{-q)\OmB{p-q)) 



ds 



xdx I da 



hK^j^ dl3hK{l-a,a-f3,f3)- ' dphxil - a,a - (3, (3)\0=o^ 
(«-m,)(l--^)-(l-^) 



K 



dl3{<t)\\ (1 - a, a - /3, /3) + (/>_L (1 - a, a - /3, /3)) 
(1 - a, a - /3, /3)|^=o + h0--a,a-(3, /3)|^=o; 



-rubfK 



d/?(/>i(l-a,a-/3,/3)- 
m%a + m^^ (1 — a) — m; 



dl3(j)±{l - a, a - /3,/3)|/3=o 



Exp{ 



{Xco{p)K{-q)\OmB{p-q)) 



Ql 



n+l 



ds 



dx I da 



^2 



1 



(46) 



1^ d(3(l>3K{l-a,a-(3,(3)- ° rf/3(/.3if(l - a,a - /3,/3)|/3=o| 



(i/3(0|| (1 - a, a - /3, /3) - 20i (1 - a, a - /?, /?)) 
d/3(</)|| (1 - a, a - /3, /3)|^=o - 2<^i(l - a, a - /?)|/3=o)| 



{x-(l-.)(l-2|)} 



mgO! + w^^^ (1 — a) — , / 



M2a 



-) 



n+l 



1 



a/3(l - a;) 



(47) 



In performing the /3 integral J^", we need only the values of the light-cone distribution amplitudes (p3K, ^||, <f>±, 
(/>l|, 4>± and their derivations at zero momentum fraction i.e. /3 = 0, there are no problems with negative partons 
(quarks and gluons) momentum fractions. The analytical continuation of to its positive value ends up with an 
unavoidable theoretical uncertainty, if only a few terms of the Taylor series are taken, smaller jad (In the Chiral 
limit m-c = 0, \s/{s — m%)\) with greater precision, however, with the re-summation to all orders of /? in Eq.(45), the 
assumption of quark-hadron duality is still applicable in the case of heavy meson final states. This procedure ensures 
the disappearance of the unphysical momentum k from the ground state contribution and enables the extraction of 
the physical matrix elements {r]c{p)K{—q)\O{0)\B{p— q)) and {xco{p)K{—q)\0{0)\B{p — q)) due to the simultaneous 
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conditions, = m\ and {p—qf' = m\. The light-cone distribution amphtudes <l>iK, <^||, 4>\\ and are analytical 
functions, as well as only numerical values are concerned, we can analytically continue into the timelike region, 
= TTig, if the unphysical negative momentum fractions of the partons (quarks and gluons) are taken for granted, 
we can then take the integral J" rf/3G(s, x, a, /3) directly, the numerical results will not make difference. In the channel 
B mr, the duality region is about = 0.7GeV^ < \P^\ ~ mf , the quantity s/{s - P^) = -s/(P^(l - s/P^)) is 
small, and the expansion in terms of Taylor series of s/P^ converges quickly, the results will not depend significantly 
on the end point values of the light-cone distribution amplitudes. As for the channels B — > rjcK, XcoK, the duality 
regions are about s^^ = llGeV"^ and s^^^ = 13GeV^, the values of s/\P^\ are about 44% and 52%, respectively. 
The expansion in terms of Taylor series of /3 for the light-cone distribution amplitudes converges more slowly and 
more terms should be taken into account, for example, \ac\ < 36% in the B rjcK channel and \ac\ < 55% in the 
B — > XcoK channel; the main drawback of this approach is the significant dependence on the end point values of those 
light-cone distribution amplitudes. 

IV. NUMERICAL RESULTS 

Next, we choose the input parameters entering the light-cone sum rules before giving numerical predictions on 
the nonfactorizable soft contributions. 

The parameters enter the decays B — > rjcif, XcO^ are taken as ms = 5.28 GeV , /b = 180 ± 30 MeV, 
rub = 4.7 ± 0.1 GeV, ss = 35 ± 2 GeV^ m^^ = 3.0 GeV, /^^ = 0.35 GeV , rric = 1.25 ± 0.05, = 11 ± 1 GeV^ 
m^^„ = 3.41 GeV, /^.o = 0.36 GeV, s^co = 13 ± IGeV^, and fx = 0.16 GeV [25]. The value of decay constant 
f^^g is taken from Ref. [26] and new estimation based on the nonperturbative approach of coupled Schwinger-Dyson 
equation and Bethe-Salpeter equation is in preparation. For the coefhcients of the twist-3 and twist-4 kaon meson 
light-cone distribution amplitudes, we can make an approximation /s^ ~ fsK, S]^ — and take fsK = 0.0026 GeV, 
S'^iHb) = 0.17 GeV, where jib = A/m| - ~ 2.4 GeV [9,23,24]. 

Here we will take a short digression to discuss the duality regions. In the B — > rjcK channel, as the axial-vector 
current = 07^750 in stead of pseudoscalar current J5 = ci'y^c is chosen to interpolate the ric meson, we must be 
careful in choosing the duality region to avoid possible pollutions from the rjc{2s) and Xci mesons with the same 
quantum numbers as the interpolating current. The masses of those two mesons are about m^^(2s) = 3.6GeV, rrixai = 
3.5GeF, and the widths of those mesons are narrow, we can choose the duality region to be s^^ = 11 ± 1 GeV^. In the 
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B — > XcqK channel, due to narrow width of the XcO meson, we can choose the duahty region to be S;^^^ = 13± IGeV"^ 
to avoid possible pollutions from the excited and continuum states. Furthermore, larger s means larger jad, more 
Taylor series of j3 have to be re-summated, heavy dependence on the end point values of the light-cone distribution 
amplitudes. In those two channels, the variation of the duality thresholds s^^ and s^^^ can lead to large uncertainties. 

The parameters n and must be carefully chosen to guarantee the excited and continuum states to be sup- 
pressed and to obtain a reliable perturbative QCD calculation. The stable region for the Borel parameter is 
found in the interval = 10 ± 2 GeV^ which is known from the B channel QCD sum rules [9]. In the charmonium 
channels, we usually perform n-th derivative and take n-th moment sum rules to satisfy the stability criteria [25] . 

For the decay B ^ r]cK, the calculation is rather stable on the range n = 3 — 7. Qo is parameterized by 
Qq = 4mg^, where ^ is usually allowed to take values from to 1 while the best interval is ^ = 0.3 — 1. In the 
following, we write down the numerical value for the nonfactorizable soft matrix element, 

{r]c{p)K{q)\d{0)\B{p + q)) = 0.035 ± O.OlOGe^^ , (48) 

the largest uncertainties come from the variations of the mass of the c quark. 

Taking into account the next-to-leading order Wilson coefficients calculated in the naive dimensional regulariza- 
tion scheme [15] for /x = m^(TOfe) = 4.40 GeV and A^^ = 225 MeV, 

CiinTbimb)) = 1.082, C2{rm{mb)) = -0.185, CsinTbimb)) = 0.014, 
C4(m^(m6)) = -0.035, C5(m6(TO6)) = 0.009 , C6(rn:^(TO6)) = -0.041 , (49) 

here we have neglected the Wilson coefficients C7, Cg, Cg, Cio in numerical calculation due to their small values, finally 
we obtain the numerical relation between the contributions from the factorizable and nonfactorizable matrix elements, 

{2KbK*.Ci(M) - ^VtbVt* [C4(m) + Cnim {f,eip)Kiq)\dmBip + q)) 



0.11 ±0.04. (50) 



From above expressions, we can see that the nonfactorizable soft contributions are considerable and they must be 
included in analyzing the branching fraction. A rough estimation shows that the theoretical branching fraction will 
increase to about 1.15—1.30 times as the naive factorization result. Although there are still large mismatches between 
the theoretical and experimental values, we can say that the theoretical prediction is considerably improved. The 
consistent and complete QCD LCSR analysis should include all the contributions from 0{as) corrections, however, 
the calculation is cumbersome and we prefer another article. 
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For the decay B —^ XcoK, the calculation is rather stable on the range n = 5 — 10, the parameter ^ is usually 
allowed to take values larger than 1 while the best interval is ^ = 1.2 — 2.0. The nonfactorizable soft hadronic matrix 

element in the B XcqK decay is 

{Xco{p)Kiq)\dmB{p + q)) = 0.22 ± O.OSGeV^ , (51) 

and the corresponding branching fraction is (1.0 ± 0.6) x 10~^ which is smaller than the experimental data (6.0 ± 
2.1) X lO"** [14]. From Eq.(48) and Eq.(51), we can see that value of the nonfactorizable soft hadronic matrix element 
{ric{p)K{q)\O{0)\B{p+q)) is about 15% of the {xco{p)K{q)\0{0)\B{p+q)) , and contributes to the branching fractions 
with magnitude about 10~^ . Take into account for the contributions from the factorizable hadronic matrix elements, 
the theoretical predicated branching fraction is about 10% of the corresponding experimental data for the decay 
B — > rjcK. In this article, we take into account only the nonfactorizable soft contributions of twist-3 and twist-4 , 
while higher twist contributions are neglected. Although there are large uncertainties due to discarding the higher 
twist contributions and varying the duality thresholds, we can estimate the order of magnitude of the nonfactorizable 
soft effects at least. 

V. CONCLUSION 

In this article, we have analyzed the contributions from the nonfactorizable soft hadronic matrix elements 
{r]c{p)K{q)\d{0)\B{p + q)) and (xco{p)K{q)\d{0)\B(j> + q)) to the decays B rj^K.XcoK with the effective weak 
Hamiltonian. As the QCD-improved factorization approach breaks down in the B to charmonia decays, the contribu- 
tions from the soft gluon exchanges will signalize themselves. Our numerical results show that their contributions are 
considerable and should not be neglected for the decay B rjcK. Although there are still large mismatches between 
the theoretical and experimental values, the theoretical predicted branching fraction is considerably improved (about 
(15 — 30)%). As for the decay B — > XcoK, there are no factorizable contributions, the nonperturbative contributions 
from the nonfactorizable soft matrix elements are about five times smaller than the experimental data. The consis- 
tent and complete QCD LCSR analysis should include all the contributions from 0{as) corrections and higher twist 
contributions. 
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